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Abstract 

Using thermal field theory, we derive simple analytic expressions for the 
spectral density of photons in degenerate QED plasmas, without assuming 
the usual non or ultra- relativistic limit. We recover the standard results in 
both cases. Although very similar in ultra-relativistic plasmas, transverse 
and longitudinal excitations behave very differently as the electron Fermi 
momentum decreases. 
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Cores of white dwarves and red giant stars are typical examples of de- 
generate QED plasmas, where the Fermi energy is much greater than the 
temperature so that the Fermi distribution function can be well approxi- 
mated as a step function. The vast majority of the literature describes the 
particle scatterings in these systems in the non-relativistic limit . But the 
electron Fermi momentum, kp, is typically of the same order as the electron 
mass. On the other hand, it is well known that the properties of light inside 
matter quite differ between the non and the ultra-relativistic limit 0, §]. 
Indeed, at low densities, longitudinal time-like excitations oscillate at fixed 
frequency u = uq, where ujq = e^Ne/rrie is the plasmon frequency 0. This is 
not so at very high densities, where the dispersion relations acquire a much 
more complicated form [Q. Also, space-like excitations seem very different 
between the two cases. 

Thermal Field Theory is a unique tool for describing a system of multi 
interacting particles in thermal equilibrium at a fixed temperature T and 
chemical potential 0, Hj. With the recent developments in the high tem- 
perature limit 0, this theory is now settled on a sound basis. Concerning 
degenerate plasmas, i.e., fi ^ T, there exists some work related with the 
emission of hypothetical particles such as axions, but dealing essentially with 
ultra-relativistic plasmas . We have been able to use the same methods in 
the general case where no approximation is made on the electron mass. The 
analytic expressions turn out to be very simple and easy to manipulate, and 
this is the main motivation for presenting our results here, separately from 
the applications we have in mind, like plasmon decay and axion emission 
i 1,011. 

The essential mathematical object we shall manipulate is the photon po- 
larization tensor, together with its finite density corrections, from which one 
can derive the photon spectral density. For illustration, we use the imaginary 
time formalism Q (our results can be derived in the real-time formalism too 
[^]). The calculations are slightly complicated by the fact that we consider 
a massive fermion However it turns out, as we show below, that the 

so-called "hard-loop" approximation can be taken in a consistent way, 
leading to simple expressions. This happens when the photon momentum 
is much smaller than the scales associated with the electron, i.e., the elec- 
tron chemical potential fi and the Fermi momentum kp (here we adopt the 

notation of a relativistic chemical potential, /i = ykp + m?). 
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The one loop contribution to the photon polarization tensor is 

d^p Tr -ff^ip - m)Y{P' - m) 



(27r)3 (P2 _ ^2) (p/2 _ ^2) 



where po = ioJn + fJ', = 7rT(2n + 1) and P' = P — Q. 

The tensorial structure is usually decomposed into transverse and longi- 



tudinal modes 11 



n 
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and, although thermal field theory is Lorentz-covariant, a preference is given 
to the plasma rest frame, where Pt and P^ are given by 
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Let us now give some details about the calculation of the transverse com- 
ponent. After performing the energy sums and taking the T = limit we 

get 
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Note that in the degenerate limit fi ^ T, the Fermi-Dirac distribution func- 
tions can be approximated as step functions which greatly simplifies the 
energy integrations. 

From now on we will ignore the /x-independent part of the polarization 
tensor (absorbed in the renormalization and subleading). Also we restrict 
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our study to the /i > case. We therefore obtain 



nT(go,g) = 
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We make our approximations at this level. In an equivalent way to the hard 
thermal loop approximation of Braaten and Pisarski, we study the case where 
the external momenta q and go are much smaller than the Fermi momentum. 
In the ultrarelativistic limit this reduces to the usual q, qo -C fi. The most 
important contribution to 11 comes from the region where p is of the order 
of kp- Then the following approximation will be valid 
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In the same limit, we can ignore the term in eq. (5) and the difference 
between 6'-functions in eq. (5) can be written as 



^(/i - Ep) - eif^ - E'p) ^ (Ep - - Ep) 

Neglecting higher orders in we obtain the final result 
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The calculation of 11^ can be performed in an identical manner. We 
obtain 



^L{qo,q) = -^fJ'kp 1 - 



lnqo , fiqo + kpq 

1 — m ; — 

2 kpq fxqo — kpq 



(9) 



As can be seen from eqs. (8) and (9), one- loop corrections start to be im- 
portant when q, qo < e^/fikp, and one should then use a resummed photon 
propagator M. 
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The dispersion relations are simply obtained from the expressions above. 
In the ultrarelativistic limit, the usual relations are recovered |10[- In 
the non-relativistic case, the dispersion relations are usually obtained from 
kinetic theory and it is interesting to see how we recover them using the 
methods of quantum field theory. Throughout this paper we have considered 
a relativistic expression for /i, that is fi = m + fip where /^i? is the standard 
Fermi energy. Therefore in the non-relativistic regime, kp <C n, and the mass 
is the highest scale we have. In a usual dispersion relation q will always be 
smaller than qo, and we can thus conclude that qkp ^ fiqo- This observation 
allows us to approximate the logarithms in eqs. (8) and (9). We find the 
following limiting values for the polarization tensor 



^dqo^q) 




^ ^ S fqkp 
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(10) 
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where we have introduced the plasmon frequency 



e2 4 e^Ne 
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which is just the relativistic generalization of the usual formula. Then, 
the dispersion relations are deduced by solving the pole equation Q'^ — 
Ren(go, q) = and we obtain 
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for the transverse oscillations and 
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for the longitudinal. These results agree perfectly with those obtained using 
kinetic theory 0. We see that the transverse oscillations have particle-like 
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dispersion relations while the frequency of longitudinal oscillations is almost 
independent of the momentum. 

In Fig. 1, we have plotted the dispersion relations for a typical value of 
the Fermi momentum encountered in the plasma core of red giant and white 
dwarf stars. As these two systems are just between the non and the ultra 
relativistic limits, the relations shown in Fig. 1 are those to be used, for 
instance, in the plasmon decay process 0. 

Whether it is more convenient to use thermal field theory or kinetic theory 
to derive these expressions for QED is probably a matter of taste. However, 
this is not the case for QCD which is a non-abelian gauge theory and where 
one is forced to use our method. As a matter of fact, the results for the 
gluon dispersion relations in a degenerate quark-gluon plasma are identical 
to those presented here, apart for some color factors. One should remember, 
however, that the quarks are mostly ultra-relativistic as the critical density 
for the transition from hadronic matter to quark-gluon plasma is much higher 
than the quark masses. 

In the calculation of emision rates of dense systems, not only the disper- 
sion relations, but also the imaginary part of the polarization tensor are of 
interest This imaginary part gives rise to Landau damping like contri- 



butions |12]. From the general expressions given above (eqs. (8) and (9)) 
we see that both the longitudinal and the transverse components develop an 
imaginary part whenever Qq is smaller than {kp/ ^)q: 

ImnT(go,g) = -i^^^f (l-(^) )^(^^^-'"5o); (15) 

Imnago,g) = ^^^''^^il-'^^^e{kpq-^^q,). (16) 
2t[ q \ q^ j 

Notice that the scales are different for the two modes. 
The spectral density for space-like excitations 

/ ^ = _1 Imn(go,g) 

/'Wo,?; 7r(Q2_Ren(go,g))2 + (Imn(go,g))2' ^ ' 

is plotted in Fig. 2 for a Fermi momentum kp = 400 keV. The static limit is 
worth emphasizing 

RenT(go = 0,g^O) = 0; (18) 
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Reni(go = 0,g-^0) = — /xfc^ = 4, (19) 



77^ 



where kn is the inverse Debye length corresponding to the well-known phe- 
nomenon of screening of static electric fields . Again, the expression given 
here is the relativistic generalization of the usual formula. Transverse exci- 
tations are not screened for any range of kp and fi. 

Despite this screening effect, one can see in the different figures presented 
here that the longitudinal spectral density has a much broader spectrum 
than the transverse one, at least in the non-relativistic limit. This originates 
from the different scales associated to the two modes. As a consequence, 
the scattering of particles will be a priori favored through the exchange of 
longitudinal photons |jl3|. Our result also illuminates the "form factor" used 
by several authors |T^ . In the language of kinetic theory, this form factor is 
related to the electric permittivity by F = (g^/fc|))(l — l/e/,) while in the 
language of thermal field theory the spectral density is pi = (l/g^)Im(l/eL) 
(in the static limit). Despite the fact that the calculations are very different in 
the two approaches, we can make some simple comparisons as the scales used 
in the problem are the same. In the static limit the photon propagator is just 
l/(?^ + ^d) cind one would naively guess that the cross section of a particle 
interacting with a charged target through the exchange of a longitudinal 
photon would behave as l/kjy However, as Imlli ~ kj-, for values of 

go — kpq/fi (which is a static point), the overall behavior is rather like l/kj^, 
which is the result found in [^. We find a justification for another point: 
in his initial work |T^, Raffelt did not take into account the damping part 



of the polarization tensor into the denominator of the propagator, obtaining 
a form factor F = q'^/{q^ + /c|)). However, it is clear from eq. (17) that if 
ImHi, ~ then it must be taken into account in the denominator. This 
point was indeed realized by Raffelt and Seckel [|T31. The obvious conclusion 



is an additionnal suppression of the scattering rates. 

In sharp contrast, these considerations do not apply in the ultra-relativistic 
limit. Magnetic and longitudinal scales become identical and both modes are 
relevant when one considers a scattering rate. Furthermore, one should not 
forget that static transverse photons are not screened 0. 

In conclusion, using thermal field theory, we have obtained simple analytic 
expressions for the polarization tensor at a finite chemical potential /i for 
all range of the Fermi momentum. The relativistic generalization of the 
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plasmon frequency is Uq = e^iVe//i and of the Debye mass kj^ = (e^/vr^)/^^^ 
(these relations can also be used when the fermion is a proton or an ion). 
In agreement with the physical picture, we have found that all collective 
effects (quasiparticle modes and Landau damping) occur when the photon 
momentum g, go is niuch smaller than the Debye mass k^- As we have also 
imposed the condition g, go kp (hard loop approximation), these effects 
are well described by our formulae only when efi <^ kp, which is the case for 
densities p > 10^ g/cm^. In any case, in a degenerate plasma, the photon 
cannot have an energy much higher than the temperature, either because of 
a Boltzmann suppression factor, or because of Pauli-blo eking effects which 
favor the soft photons exchange [Q. As yU, fci? 3> T, one can see that the 
expressions presented here (eqs. (8) and (9)) are rather general and will be 
of practical use, especially in stellar systems where neither the non nor the 
ultra relativistic approximation is accurate. 
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Figure Caption 



Fig. 1 Dispersion relations for timc-likc excitations for a Fermi momentum 
kp — 400 keV. The solid line represents the transverse mode and the 
dashed line the longitudinal mode. 

Fig. 2a The spectral density for transverse space-hke excitations for a Fermi 
momentum kp — 400 keV. 

Fig. 2b The spectral density for longitudinal space-hke excitations for a 
Fermi momentum kp — 400 keV. 
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